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Abstract 

For assemblies of superconductor nanograms, the magnetic response is analyzed as a function 
of both temperature and magnetic field. In order to describe the interaction energy of electron 
pairs for a huge number of many-particle states, involved in calculations, we develop a simple 
approximation, based on the Richardson solution for the reduced BCS Hamiltonian and applicable 
over a wide range of the grain sizes and interaction strengths at arbitrary distributions of single- 
electron energy levels in a grain. Our study is focused upon ultra-small grains, where both the 
mean value of the nearest-neighbor spacing of single-electron energy levels in a grain and variations 
of this spacing from grain to grain significantly exceed the superconducting gap in bulk samples 
of the same material. For these ultra-small superconductor grains, the overall profiles of the 
magnetic susceptibility as a function of magnetic field and temperature are demonstrated to be 
qualitatively different from those for normal grains. We show that the analyzed signatures of 
pairing correlations are sufficiently stable with respect to variations of the average value of the 
grain size and its dispersion over an assembly of nanograins. The presence of these signatures does 
not depend on a particular choice of statistics, obeyed by single-electron energy levels in grains. 

PACS numbers: 74.78.Na,74.25.Ha,74.20.Fg 
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I. INTRODUCTION 



Recent intensive theoretical studies of pairing correlations in small superconductor grains 
(see Ref. Q for a review) are inspired, to a great extent, by the experiments of Ralph, Black 
and Tinkham 0, 0, 1^ , who measured tunneling spectra of individual nanosize Al grains. For 
grains with an even number of electrons and the mean value of the nearest-neighbor spacing 



for single-electron energy leye^ 
Ab, the measurements u 



Sj d, significantly smaller than the bulk superconducting gap 
revealed a gap between the ground state and the lowest 
excited state of a grain. The width of this spectroscopic gap is larger than or comparable 
with 2Ai, reflecting the fact that the energy difference between the fully paired ground state 
and the lowest excited state, which involves one broken pair, is at least as large as the 
energy cost for pair breaking. The lowest excitation energy contains two contributions: one 
contribution (~ d) is given by the energy spacing between the corresponding single-electron 
energy levels and just coincides with the lowest excitation energy for a normal grain, while 
another (~ 2Af,) corresponds to the change of the pairing-interaction energy due to pair 
breaking. When the mean spacing d is much smaller than A;,, the lowest excitation energy 
considerably exceeds d so that a specroscopic gap in the tunneling spectrum can be clearly 
observed and unambiguously attributed to the presence of pairing correlations in a grain. 
However, in the case of ultra-small grains with d > A^, the lowest excitation energy is mainly 
determined by the energy difference between single-electron energy levels and hence pairing 
correlations can hardly be detected in a convincing way through the tunneling spectra of 
individual grains. 

For ultra-small grains with d > A^, traces of superconducting correlations in thermody- 
namic properties (specific heat, magnetic susceptibility) have become a subject of theoretical 
investigations 

mm 

. In Refs. p, laLD], pairing effects in the spin magnetization are 
studied for a particular case of grains with uniformly spaced single-electron energy levels. 
Signatures of pairing in the magnetic response of those grains were shown to persist even for 
d ^ Afe. In particular, the study of the magnetic susceptibility in the limit of zero mag- 
netic field revealed that the susceptibility of an odd superconductor grain (for short, we will 
use the terms 'even grain' or 'odd grain' for a grain with an even or odd number of electrons, 
respectively) should have a minimum as a function of temperature T at k^T ~ d, where k-Q 
is the Boltzmann constant. This minimum originates from a superposition of the Curie-like 
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contribution of a single electron, which is unpaired at T = 0, and contributions, which are re- 
lated to electrons paired at T = and increase with temperature due to pair breaking. Such 
a non-monotonous behavior is visible in the calculated magnetic susceptibility jsj of grains 
with the ratio dlAi, as large as 50. Calculations of the zero-field magnetic susceptibility have 
been extended j8| to assemblies of superconductor grains, assuming that single-electron en- 
ergy levels in grains follow the Gaussian orthogonal ensemble (GOE) distribution The 
zero-field magnetic susceptibility of assemblies, which contain only odd grains or approxi- 
mately equal fractions of even and odd grains, as a function of temperature is shown to 
manifest a non-monotonous behavior similar to that described above. The zero-temperature 
magnetic response of superconductor nanograins with the GOE statistics of single-electron 
energy levels has been investigated |9] using a perturbative approach to calculate the inter- 
action energy from the Richardson equations 11], which provide an exact solution for the 
reduced BCS Hamiltonian. While for normal grains the zero-temperature magnetic suscep- 
tibility monotonously increases with magnetic field H, a maximum at ~ d (/iB is the 
Bohr magneton) was found in the susceptibility of superconductor grains as a function of 
H. For ultra-small superconductor grains with the dimensionless interaction constant as 
large as 0.28, the calculated height of this maximum jS] is by two percent larger than the 
saturation value of the magnetic susceptibility for both superconductor and normal grains 
at if ^ oo. 

This paper is aimed to study the magnetic response of superconductor grains as a func- 
tion of both temperature and magnetic field. Our study is focused upon ultra-small grains, 
where both the mean value of the nearest-neighbor spacing of single-electron energy levels in 
a grain and variations of this spacing from grain to grain exceed the superconducting gap Af, 
in bulk samples of the same material. For the ultra-small grains, we analyze observable man- 
ifestations of pairing correlations in the behavior of the magnetic susceptibility as a function 
of T and H . Since magnetization measurements on single nanograins look problematic, we 
perform calculations for assemblies of grains. It is worth noting that an appropriate choice 
of statistics to describe the energy-level distribution depends on concrete nanostructures. 
Thus, ultra-small superconduct or g rains, p repared by ion implantation, are known to pos- 
sess a regular geometric shape [ij. Il^ ll^ . implying that single-electron energy levels can 
be degenerate and hence that they do not obey the GOE statistics, often used for supercon- 
ductor nanograins in the literature. Below, we consider two cases: (i) parallelepiped-shaped 
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grains [15|, where single-electron energy levels obey the Poisson statistics and (ii) grains 
with the GOE statistics of single-electron energy levels. In order to make the analysis of 
the magnetic response possible for large assemblies of nanograins at non-zero temperatures, 
when a large number of many-particle states are involved in calculations, we develop an 
approximation for the interaction energy of pairs in a nanograin. This approximation, de- 
rived on the basis of the Richardson equations combines the advantage of simplicity 
with the possibility to analyze the effects related to an extremely high sensitivity of the 
interaction energy to a specific distribution of single-electron energy levels. 

The rest of this paper is organized as follows. In Section II we describe the model used 
to calculate the magnetic response of superconductor nanograins. In Subsection II. A, we 
consider the Richardson solution for the reduced BCS equation as applied to the energy of 
pairs in the absence of a magnetic field. Energy levels in nanograins in the presence of a 
magnetic field and the spin magnetization are analyzed in Subsection II. B. In Subsection II. C 
we develop an approximation for calculating the interaction energy of pairs in a nanograin. 
In Section III the proposed approximation is used to study the magnetic susceptibility of 
superconductor nanograins with different parity of the number of electrons as a function of 
magnetic field and temperature. We examine also the effect of a particular choice of the 
single-electron energy-level statistics on the magnetic response of superconductor nanograins. 
Further, we consider a dependence of the magnetic susceptibility on the interaction strength 
as well as on the grain size and its variations within an assembly of grains. Conclusions are 
formulated in Section IV. In Appendix A we analyze - as an instructive example - a grain, 
where in the absence of a magnetic field single-electron states in the interaction band have 
equal energy. In Appendix B the accuracy of the proposed approximation for the interaction 
energy is checked by comparing the magnetic susceptibility, calculated on the basis of this 
approximation, with the numerically exact results for some particular cases. 
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II. FORMALISM 



13) 



A. Energy levels at zero magnetic field 

For a nanosize grain, the BCS pairing Hamiltonian can be written as (see, e. g., Refs. Q, 



H e^al^tti^ - XdY^ al^al^aj^aj^, (1) 

where the operator i^ja) creates (annihilates) an electron with the spin a in the i-th 
spin-degenerate single-electron state. We will assume that the index i {i = 1,2,...) labels 
these states in the order of increasing energy Si. The second term in Eq. is the interaction 
Hamiltonian. The sum in this term is over the set X of / states inside the energy interval 
{e-p — /k^D, + ^t^o), which will be referred to as the interaction band, with Tiwd, the Debye 
energy, and ep, the Fermi energy. The interaction strength is a product of the mean energy- 
level spacing within the interaction band, d = 2hujD/I, and the dimensionless interaction 
parameter A. 

The electrons occupying levels outside the interaction band are straightforwardly de- 
scribed by the first term in the r.h.s. of Eq. (^. Therefore, we will focus now on the 
electrons, which reside in the interaction band. Let the interaction band be populated by 
2n + b electrons, of which b electrons are on singly occupied levels (the set of these levels, 
blocked to the pair scattering, is denoted as B). The electrons on singly occupied levels do 
not participate in the pair scattering and their contribution to the energy of the iV-electron 
system under consideration is simply XligB^j- The remaining 2n electrons form n pairs, 
which are distributed over the set W = I\ BoiI — b unblocked levels. Richardson 
showed that the energy of these pairs is given by the sum of n parameters Ei, . . . , En, which 
are non-degenerate roots {E^ ^ Ey for all v) oi the following set of n coupled equations: 

— -Y ^ + y ^ = 0, z/ = l,...,n. (2) 

If the interaction parameter A reduces to zero, each solution of Eq. (0) smoothly evolves 



into a certain set of bare pair energies In other words, for a given set U of unblocked 
bare energy levels in the interaction band, a specific choice of a subset M oin levels, doubly 
occupied in the limit A ^ 0, uniquely determines a solution of Eq. (j21) at A 7^ 0. Due to such 



a one-to-one correspondence between many-electron states at A = and those at A 7^ 0, the 
states in a superconductor grain can be classified according to the hmiting distributions of 
electrons over bare energy levels for A — 0. 

Let the total number of electrons in a grain be 2K + P, where the non-negative integer 
K coincides with the number of doubly occupied bare energy levels in the ground state at 
A = in the absence of a magnetic field, while the parity parameter P takes the values 
or 1 for, respectively, an even or odd grain. For the energy levels of this {2K + P)-electron 
system in the absence of a magnetic field we use the notation Eq^, where K — k gives the 
number of doubly occupied bare energy levels at A = 0, while q labels a specific distribution 
of K — k doubly occupied levels and {2k + P) singly occupied levels within the set 
Let T'('?'^) and S^'^'^^ denote, respectively, the sets of doubly and singly occupied energy levels 
Ei for the state with given q and k. The parameters q and k uniquely determine also the 
sets U^''''\ and J\f^i''^ of singly occupied (blocked), unblocked, and doubly occupied (at 
A = 0) bare energy levels in the interaction band, as well as the number of levels in these 
sets: b^'^''\ u'^'^^^ = I — h^'i^\ and n'^'^^\ respectively. For fixed values of q and fc, only one 
solution of equations (121) reduces to the set of bare energy levels {2e.(qk)} (/i = 1, . . . ,n*^'^'^-*), 
where ifl^^ are the elements of the set M'^'^^\ The roots, which belong to this solution, can 
be denoted as E^^''^ (/i = 1, . . . , n^'^^^). In parallel to the set of e'^^''^ (/i = 1, . . . , n^'*'^)), it is 
convenient to introduce the equipotent set of parameters e\''^'^ {i G A/'^'^'^-'), determined by 
the relation E%% = Ej^^'"'^ (/x = 1, . . . , n^^''^). In the set {E^"'"^}, each root of equations 
is labelled in accordance with its limiting value at A — > 0, that is 2ei (i G A/''-'^'^-*). So, the 
energy Egk takes the form 

where the interaction energy for n^''^^ electron pairs, which reside in the interaction band, is 

with E^'^''\ the roots of the equations 

Tl-T ^+ y -7T7^-7TY = 0' t e X^^'K (5) 
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Let us label with q = the lowest (2i^ + P)-electron state with a given number k. In this 
state at A ^ 0, 2{K — k) electrons fully occupy the lowest {K — k) energy levels e^, while 
the next {2k + P) energy levels Si are singly occupied by the remaining {2k + P) electrons. 
The corresponding energy is 

K-k K+k+P 
i=l i=K-k+l 

It is obvious that the interaction energy i?Q^ differs from zero only for k < 72^°°^ , where n^^'^^ 
is the number of electron pairs in the ground state, which has the energy 

K 

Eoo = 2Y,ei + PeK+i-El,. (7) 

i=l 

From Eq. 0, one can easily derive the relation 

Afe-Eofc = ^^kEl^ + {ex+k+p+i — ^K-k) , (8) 

where the forward difference of any function of an integer variable k is defined as A^Fk = 
Fk+i — Fk- It is worth mentioning that for k < n'-'^^-' the quantity {—A^Eqi^), which has the 
meaning of the pair-breaking energy, is strictly positive. Therefore, as follows from Eq. ©, 
the energy Eok is a strictly increasing function of k for k < n^^^\ (At k > n^^^^ the forward 
difference A.kE\^ may take zero value if the energy levels Si are degenerate.) 



B. Spin magnetization 

When analyzing the magnetic response of nanograins, we assume these grains to be flat 
and subjected to a magnetic field parallel to the grain base. For an in-plane magnetic field, 
the orbital magnetization of electrons in a fiat nanograin can be neglected Q,!^; 21 1. In a 
magnetic field H, the energy of a {2K + P)-electron state depends on the spin orientation 
of electrons on singly occupied energy levels Si. 

Eqkl = Eqk -{2k + P- 21) fiBH. (9) 

Here / (/ = 0, ...,2k + P) is the number of spin-down electrons on singly occupied energy 
levels Ei (we use the terms 'spin-up' and 'spin-down' for electrons with spins, which are, 
respectively, parallel and antiparallel to the applied magnetic field), and the Lande factor 



for the electron spin is assumed to equal 2. The spin magnetization of electrons in a state 
with the energy (jH)), does not depend on q: 

Mki = {2k + P - 21) fiB. (10) 
At a temperature T, the spin magnetization of the grain can be expressed as 

= I E E cr^'^f- E (-ff ) . (11) 

fc=0 1=0 9 ^ ^ ^ 

with the partition function 

K 2k+P / p \ 

2 = EEc^?'""E-p(-tf)' (12) 

The presence of the binomial coefficients Cf^^^ in Eqs. (fTT|) and (fT^ takes into account 
the degeneracy of the energy level Egki with respect to permutations between spin-up and 
spin-down electrons on singly occupied energy levels Ei. 

To facilitate the further analysis, let us first consider magnetization of a nanograin at 
zero temperature. At T = 0, magnetization is determined by the polarization of electron 
spins in the ground state. To find the ground state at a given magnetic field, one has to 
minimize the energy Eqki with respect to the numbers q, k, and /. As obvious from Eq. Q, 
in the ground state the number of spin-down electrons on singly occupied levels must be 
zero: / = 0. We also note that only the first term in the r.h.s. of Eq. Q does depend on 
q. At a fixed number k, the minimal value of this term is Eqi. [see Eq. ©]. Therefore, the 
magnetization of the grain at T = is given by the value M^o with the number k, which 
minimizes the energy 

Eoko = Eok - {2k + P) fiBH . (13) 

Since the energy Eok has been shown to be a non- decreasing function of k (a strictly in- 
creasing function of k for k < n^'^^^), at if the number of singly occupied energy levels 
Ei in the ground state goes to P and, correspondingly, the magnetization is given by the 
expression Mqo = /xb-P- However, as seen from Eq. ()13|) with rising H an increase of k can 
become energetically favorable. Each jump of the value km, which minimizes the energy 
Eoko {H), results in an abrupt increase of the magnetization. 
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Under the condition A^Eoko > or, equivalently, 

^lEok > 0, (14) 

the energy Eoko (H) as a function of k has only one minimum on the interval [0, K], while 
the value km is the minimum integer k, which obeys the inequality AfcE'ofco > 0. In this 
case, an increase of H results in consecutive jumps of the value k^, from A; to A; + 1, where 
k = 0,1,2,.... The magnetic field Hk, which corresponds to the change of km from k to 
k + 1 is determined by the equation A^i^ofco = 0. From Eq. |T^. we find 

AfcE'ofco = AkEok — 2fXBH (15) 

and, using Eq. (jH)), 

Hk = - — (— Afc-Eofc + ^K+k+p+i — £K-k) ■ (16) 

The magnetic field Hk given by Eq. (fT6|l is a non-decreasing function of k, provided that 
the inequality (fT^ is obeyed. From Eq. (jS)) we obtain the equation 

A-lEok = —AkEli^ + {SK-k — ^K-k-l) + i£K+k+P+2 — ^K+k+P+l) , (17) 

which clearly shows that inequality p4j] is always satisfied for normal grains, where the 
interaction energy is identically zero. However, for superconductor grains, the sign of A^E^k 
is not obvious. On the one hand, for k < n^'^'^\ an increase of k reduces the number of 
electron pairs n^^^'> and, simultaneously, enlarges the number b^^''^ of singly occupied energy 
levels Si in the interaction band, which are blocked to pair scattering. Both a decrease of n^^*^^ 
and an increase of b^^''^ weakens the pairing interaction. For this reason, the pair-breaking 
energy (— A^i^Qj^.) decreases with increasing k and, correspondingly, the values of (— A|i?Q^) 
are negative. On the other hand, for sufficiently small values of the energy differences 
i^K-k - £K-k-i) and {eK+k+p+2 - Sx+k+p+i), the sign of AlEok in Eq. ((T7j) is just the same 
as the sign of (— A^i?g^). Thus, depending on a specific energy spectrum {si} as well as on 
particular values of A and k, the inequality (jl4j) may be violated in superconductor grains 
(an instructive example, which corresponds to the case of A^EQk < 0, is given in Appendix 
A). A violation of the condition (jl4j) means that the energy E^ko as a function of k can 
have two or more minima on the interval [0,7^]. Since the relative depth of those minima 
depends on the magnetic field, with rising H an abrupt increase of km by a value Ak > 1 
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is possible. As clear from Eq. ()15p and the definition of Hk-, a given value k can correspond 
to the ground state only for H < Hk- Therefore, an abrupt change of km from km = k to 
may happen only at a magnetic field Afc, which satisfies the inequality 

Hk,Ak < Hk- (18) 

As follows from the above considerations, at T = the magnetization M of a single 
nanograin forms a staircase, each jump of M resulting in a (5-like peak of the magnetic 
susceptibility x- From Eqs. (fT^ and (fTSj) . it is obvious that positions of these peaks are 
strongly affected by a specific shape-dependent configuration of the energy levels in the grain. 
Since this configuration is hardly controlled in experiment, a possibility to convincingly 
reveal pairing correlations through the susceptibility measurements for a single nanograin 
looks unlikely (apart from obvious difficulties of those measurements as such). Therefore, 
we analyze below the magnetic susceptibility for assemblies of grains, where some general 
trends related to the pairing interaction can be manifested. 

For an assembly of grains, we define the magnetic susceptibility as the forward difference 
quotient 

where angular brackets (...) stand for the average over the assembly of grains. The above def- 
inition of the magnetic susceptibility is more convenient than the derivative (x) = d{M)/dH 
because at T ^ the latter becomes a set of 5-like peaks for any finite number of grains. 
In the calculations below, the values AH ~ 10^^{d)/^B are used. We introduce also the 
normalized magnetic susceptibility 

C. Approximation for the interaction energy 

An exact solution of Eqs. (0) can be obtained only by numerical computation, which 
becomes more and more time-consuming with increasing the number of single-particle energy 
levels and the number of electrons in the interaction band. At the same time, calculations 
of the spin magnetization for a nanograin require to compute interaction energies for many- 
electron states, whose number rapidly increases with increasing temperature. Under those 
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circumstances, an exact calculation of the spin magnetization seems problematic, especially 
for large assemblies of grains. Because of that, here we introduce an approximation, which 
combines the advantage of simplicity with the possibility to catch the effects ^| related to 
the extremely high sensitivity of the interaction energy to the distribution of single-electron 
energy levels over the interaction band. 

In this subsection, we focus our analysis on the system of the electron pairs distributed 
over a fixed set U of u unblocked spin-degenerate energy levels Si in the interaction band. 
So, when discussing, e. g., the ground state (excited states) of the system, we mean the 
lowest state of these electron pairs (their collective excitations) irrespective of the state of 
other electrons, contained in the grain. For brevity, we do not explicitly indicate the indices 
q and k, which label many-electron states of the whole grain. Unless otherwise stated, we 
assume to deal with any of those states. 

We will construct our approximation for the interaction energy as an "interpolation" 
between the results for the limiting cases of weak and strong pairing interaction. 



1. Weak interaction 



The case of weak interaction is defined by the inequality 



\d <^ \ei — eA, for i,j&U, i ^ j. 



(21) 



According to Eq. (0), we have 



2ei - E, = Xd 



Xd 



2Xd 



Oc . — E- ^-^ E - — E- 



(22) 



Under condition (PT|). all the roots Ej (j G A/") in the r. h. s. of Eq. (j22I) can be approxi- 
mated by their limiting values at A — 0. Then the interaction energy becomes 



(23) 



where 



A,; = A 



1 - Xd 



(24) 
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Formally, equation ()23p looks as a sum of interaction energies for n independent subsys- 
tems, where each subsystem consists of a single spin-degenerate energy level with a single 
electron pair on it, the pairing interaction for the i-th subsystem being characterized by the 
interaction parameter Aj. The dependence of the interaction energy on a specific choice of 
a many-electron state is straightforwardly taken into account through the dependence of Aj 
on the sets U and A/" specific for each many-electron state. 

2. Strong interaction 

The case of strong interaction is defined by the inequality 

Xd^5 (25) 

with 6 = 2ma.Xi j(zi( \ei — ej\. Inequality (j^^ can be satisfied even at A < 1 provided that the 
energy levels ej in the interaction band are (quasi-) degenerate. Under the condition ((23), 
the roots Ei {i G A/") as a function of Xd demonstrate two qualitatively different types of 
behavior. While for the ground state all these roots are characterized by \2ei — Ei\ ~ Ad 
and tend to infinity at Ad ^ oo, for excited states a subset Q d M oi g roots remain finite 
at Xd oo: \2ej — Ei\ < 5 for i & Q and j G U. The roots, which belong to the subset 
are shown to represent bosonic pair-hole excitations of the system 22|. The algorithm. 



developed in Ref. |2^ allows one to count the number of those roots depending on a specific 
choice of a subset M within the set U of unblocked energy levels in the interaction band. 

Taking into account the inequalities 2\eu — Sj\ < 5 {j E U) and \2eu — Ei\ < 5 {i E Q) 
with Eu = u^^ Ylij&A^j^ '^^ convenient to rewrite Eqs. (0) for i G A/" \ ^ in the form 



OO 

1 u — 2q \ ^ 2 2mi x ^ 2mr, — Ur, . r i ^ / n 

VI ^ + V^ ^ + ^ + V " =0, ieJ\f\g, 26 

Xd 2eu-E, f^E,-E, {2eu - E,f {2eu - E,)p 



where 



/i, = 5^(2£^,-2£,)^ (27) 



J2i2s-u-E,r. (28) 
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Using Eq. ()26|). we obtain for the interaction energy the expression 
= Xd 

with ej'/ = n^^ Y,jeJ\f^j ^^"^ 



[n - g){u -n-g+1) — + — 



2m„ 



Xd 



p=2 



{xdy 



(29) 



jeN\g 



Xd 



2eu — Ei 



(30) 



A set of equations for calculating the coefficients Sp is derived by summing up Eqs. (j^Bjl . 
multiplied by {2eu — Ei)~^ (r = 0, 1, 2, . . .), over i G Af. Thus, omitting terms ~ 5"^ /{XdY 
and ~ 6/{Xd) in the sums with r = and r = 1, respectively, we obtain 



2s2mi 

n- g-[u- 2g)si + = 0, 

si — — (m — 2(7 — l)s2 = 0. 
With si and S2 determined by Eqs. and (jH^ . the interaction energy becomes 



(31) 
(32) 



= Arf (n-5()(M-?2-5( + l) 

(n - 5()(m - 5- - n)(/i2 - 2m2 



2n{eu - _^ ~ 2n)mi 



Xd 

2ml) 



{u - 2g)Xd 



+ ... 



(u-2gnu-2g-l){Xdy ^^^^ 
where the omitted terms are of the order of ~ 6'^/ {Xd)^ as compared to the first term in 
square brackets. Equation (jH^ extends the results for the energy of n electron pairs on a 
single M-fold degenerate energy level in the interaction band (see, e. g., Refs. I26ll27ll to the 
case when this level is quasi-degenerate. It is noteworthy that in the important particular 
case of half filling {u = 2n) a dependence of the interaction energy on the values Ei with 
i & Q (which in general cannot be found analytically) appears only in the terms of the order 
of ~ 5^ / {Xdy [p > 2) as compared to the leading term. 



3. Arbitrary interaction 

When deriving an approximation applicable for arbitrary values of the ratios Xd/{ei — Sj) 
G W), we start with the ground state of n pairs, when the set Af corresponds to the 
lowest energy levels from the set U. For each i E Af we define a subset Ui <ZU oi Ui energy 
levels Sj in accordance with the conditions 

— ^il < c for j G Ui, (34) 
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I^i — ej\> c for i ^ Ui. 



(35) 



The way we choose the value of the "discrimination parameter" c will be specified later [see 
the text after Eq. Let Ni denote the subset of rii doubly occupied levels within the 

set Ui. Ni = N \~\Ui. Then equations (0) can be rewritten as 



1 

Ad 



2e,; — Ei 



Ui - 



^ 2ei - 2e 



Ei 



( 



+ 



- E E 



2ei — Ei 
0, ■ 



rii 



\ ^ Asi — Ei — Ei 

1 + 2^ J 

J 



Ej — Ei 



eAf. 



(36) 



When raising Eq. to an argument for applicability of a perturbative approach similar 
to that considered in Subsection II.C.l, the differences 2ej — Ei and Ej — Ei {j ^ Ui), which 
enter the last two terms in the 1. h. s. of Eq. (|36|) . can be approximated by their limiting 
values at A — s> 0. Summing up Eqs. (fHHjl . multiplied by 25 j — Ei, over i G Af, we find for the 
ground-state interaction energy: 



El 



g.s. 



^ Xid {ui -rii + l) + ^Xid^ 

(Ej + Ei){\i - \j) - 4(Ai5i - XjEj) 



+ 9 E 



Ej — Ei 



(37) 



where 



A,: = A 



^ 2(5,- — Ei) 



(38) 



When treating Eq. (j3l|l as a criterion of (quasi-) degeneracy of single-electron states, which 
belong to the set Ui, the inequalities 2|5j — 5j| <^ |25j — Ei\ and 2|5j — 5,;| ^ \Ej — Ei\ 
are satisfied for %,i G Ui [cp. the analysis of the ground-state interaction energy under 
condition (^3)) in Subsection II. C. 2]. Therefore, the last two terms in the r. h. s. of Eq. ()37p 
can be neglected as compared to the first term. Thus, we obtain the expression 

Xidui {ui - + 1) 



El 



g.s. 



E- 



Ui 



(39) 



Expression ()39p corresponds to the approximation of (non-degenerate) energy levels Sj 
{j G Ui) by a single degenerate level. As implied from Eq. taken at g = 0, such an 
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approximation overestimates the ground-state interaction energy. Obviously, the extent of 
this overestimation increases with increasing c. On the other hand, an underestimate of the 
parameter c resuhs in oversized values of Aj, determined by the perturbative expression (|38|). 
As seen from Eq. (j38|) . these values have a tendency to diverge for c ~ Xd. For the above 
reasons, equation (j39|) is believed to provide an upper bound for the ground-state interaction 
energy. Therefore, it is natural to choose the parameter c by minimizing the value M ^ given 



The ground-state interaction energy, described by Eq. (jH^ . can be considered as a sum 
of the interaction energies for n mutually independent systems weighted with the factors 
n~^. Each of these systems contains rii electron pairs, distributed over a set Ui of degenerate 
single-electron states, and is characterized by the interaction parameter Aj. Assuming that 
the above structure of the expression for the interaction energy keeps for all the n-pair 
states, we extend our approach to the case of collective excitations of electron pairs by 
approximating the interaction energy with the expression 



In Eq. pn|l . the number of elementary collective excitations Qi is calculated in accordance 
with the algorithm applied separately to each set Ui with the subset A/^ of doubly 
occupied levels. In the limit A ^ oo, when the inequality is satisfied, the sets Ui (A/^) 
for all i E Af coincide with U (Af) so that the r. h. s. of Eq. (jlU)) reproduces the leading 
term of the result given by Eq. (jHHj) . In the opposite case described by Eq. (jUJ, both the 
sets Ui and Afi reduce to a single element, i, and hence equation (jlUjl transforms into the 
perturbative result (j23|) . 

In order to assess the accuracy of the proposed approximation, we compare in Fig. the 
exact and approximate results for the condensation energy. 



in the ground state of a nanograin. The last term in Eq. (j4ip describes the interaction energy 
for the uncorrelated Fermi ground state. The calculations are performed for parallelepiped- 
shaped nanograins with a fixed even number of electrons, 2K = 4000, and various shape, 
described by the aspect ratio : ly h = 1 '■ (1 + ^) '■ (1 + 2?7). The absolute values of k 
{i = x,y,z) are determined by the equation 2K = {lx){iy){h)ne, where rig is the electron 



by Eq. 




(40) 




E — nXd 



(41) 
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density. The single-electron energy spectrum in a hard- wall parallelepiped with sizes Ix, ly 
and Iz is given by the expression 

''^^ [l^ i^j ' Mi),nyit),nzii) = 1,2,3, (42) 

where m is the effective mass of an electron. The correspondence between i and nx{i), 
ny{i), riM) is determined by the requirement Si < Sj+i, which we impose. Henceforth, the 
values [l8| of m and ng for Al are used: m = lArrie and = 181 nm~^. In the ultra-small 
grains under consideration, the spacing d varies from 1.5 meV to 3.4 meV depending on 
the parameter rj. As expected, our approximation is seen to overestimate the condensation 
energy and, correspondingly, the interaction energy. However, the difference between the 
exact and approximate values appears to be significantly smaller than their variations with 
changing the shape of nanograins. Due to an extremely high sensitivity of energies Ei to 
the grain shape, even within a set of grains with a relatively weak dispersion of shape 
{i] G [0, 0.06]), one finds a vast variety of different single-electron energy spectra: from those 
with a single multiply degenerate energy level in the interaction band (at rj = 0) to the 
spectra with quasi-uniform distributions of levels over the interaction band. As implied 
by Fig. HJ the developed approximation provides an adequate description of the interaction 
energy in all these cases. 

While the results shown in Fig. prelate to ultra-small grains, our next goal is to estimate 
the accuracy of Eq. ()4()j) in the opposite limit of n — oo. We perform such an estimate for 
the ground state of n electron pairs in the half-filled {u = 2n) interaction band formed by 
the equidistant single-electron energy levels (^j+i — Si = d with d = huY)/n, i G XY The 
exact result for the ground-state interaction energy of such a system at n ^ cxo is [23[ 



=nLUj:>n\\ll + 



sinh 



, A, A > 1, 

hujDn X < (43) 
2exp(-2/A), A<1. 

In the case of A — *^ oo, our equation (j40p becomes exact resulting in the same value of the 
ground-state interaction energy as that given by the corresponding expression of Eq. (|43|) . 
In the case of A <^ 1, under the reasonable assumption c ^ we obtain from Eq. ()4()|1 the 
expression E^^ ^ c^nX{2}TWj))^^[l — A ln(/za;D/c)]~^, which must be minimized with respect 
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to c. The minimizing value of c appears to be c = ^Dexp(— 1/A + 1/2), which is indeed 
much smaller than hujD- Consequently, the interaction energy takes the value 



which differs from the corresponding result of Eq. by the factor e/2 ^ 1.36. 

The accuracy of the developed approximation is amazing when taking into account the 
simplicity of this approximation and the possibility to use it over a wide range of A and n at 
arbitrary distributions of single-electron energy levels in the interaction band. The accuracy 
of the results for the magnetic susceptibility, obtained on the basis of our approximation, is 
analyzed in Appendix B. 

III. MAGNETIC SUSCEPTIBILITY FOR ASSEMBLIES OF GRAINS 
A. Pairing and parity effects in the magnetic susceptibility 

We start with the magnetic response of parallelepiped-shaped grains similar to those 
considered in Subsection II. C. 3. Within an assembly of grains, grain sizes /j {i = x,y,z) 
are assumed to vary as k = + ri6), where are random numbers (—1 < < 1) 

and 6 is the amplitude of the relative size deviations. In the further calculations, we take 
(Ix) '■ {ly) '■ {Iz) = 1:7:7 and 5 = 0.025. The absolute values of k {i = x,y,z) are 
determined by a choice of the parameter Kq in the equation 2Ko + P = {Ix) {ly) {lz)ne, which 
gives the number of electrons in a grain with volume {lx){ly){h)- 

In Fig. 121 we show the normalized magnetic susceptibility xi^^ T) calculated for two 
assemblies of even grains. FiguresEfa) andEfb) display the results for normal grains (A = 0) 
and for superconductor grains (A = 0.22), respectively. In both cases, the parameter Kq is 
2000 and the corresponding average nearest-neighbor level spacing in the interaction band 
is (d) = 3.1 meV. 

As shown above, at T = 0, non-zero values of the magnetic susceptibility of a single 
normal grain correspond to jumps of magnetization, which occur at the magnetic fields 
given by Eq. (fTI)|) . For an assembly of grains at T = 0, the magnetic susceptibility, defined 
by Eq. (fTIH) . is determined by contributions of grains with Hk G {H, H + AH), where k takes 
any allowed value {0 < k < K — 1). When taking the limit of an infinite ensemble of grains 
and AH 0, the magnetic susceptibility {x{H, 0)) is proportional to the probability density 




(44) 
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of finding Hk = H (with any allowed value of k) within this ensemble. As we mentioned in 



Ref. |15|, at if = parallelepiped-shaped grains are characterized by the Poisson distribution 
of the spacing between single-electron energy levels. This means that the probability density 

Wm{s) for a particular value s of the spacing between the levels ej+m+i ("^ = 0, 1,2, ) 

and Ej in a grain with the mean nearest-neighbor level spacing d takes the value 

<°'^^°°(^) = Afll^expf-^). (45) 



ml \dJ \ d/ 

Taking into consideration Eq. (jl6|) and assuming K ^ 1, the magnetic susceptibility x for 
an (infinite) ensemble of normal grains can be approximated by the expression 

/I 2 °° 

ix) = -^J2'^2kM'^fiBH), (46) 

k=0 

where the sum over m determines the total probability of finding a spacing 2^^H between 
two single-particle energy levels, separated from each other by any even (for P = 0) or any 



odd (for P = 1) number of other levels (cp. Refs. l24l. l25l). For an ensemble of even grains 
with the Poisson distribution of the energy-level spacing, one obtains from Eqs. ()45|1 and 

ini): 



ix) 



Poisson 

" d 



1 + exp 



d 



(47) 



As seen from Fig. Efa), the results of our numerical calculations for normal grains at T — 
basically are in agreement with Eq. (j37j). Deviations of the calculated values x{H,0) from 
the behavior prescribed by Eq. are caused by the following two reasons: (i) fluctuations, 
related to seemingly random distributions of energy levels in individual grains [l^, are 
not completely suppressed even for an assembly of as many as 60000 grains; (ii) for our 
assemblies of parallelepiped-shaped grains, the value d is somewhat different for different 
grains. 

From Fig.|21^a), the magnetic susceptibility x{H, 0) for even normal grains is seen to have 
a maximum at if ^ 0. This maximum is due to the fact that the function WQ°^^^°^{2fiBH) , 
which describes the probability density for the nearest-neighbor spacing to be 2^^H , takes 
its maximum value at if ^ 0. As distinct from normal grains, in even superconductor 
grains a transition from the nonmagnetic state with = (when all the electrons are on 
doubly occupied energy levels Si) to a magnetic state with k>l always requires a non-zero 
energy for pair breaking. Therefore, at T ^ the low- field magnetic susceptibility for even 
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superconductor grains is zero. A non-zero magnetic response appears only at relatively high 
magnetic fields, when the magnetic energy becomes close to the spectroscopic gap in 

grains. 

In Fig. |2tb), one finds a well-pronounced maximum of x ^ function of H at low 
temperatures. In what follows we explain the nature of this maximum. 

First, let us consider grains, where increasing H leads to consecutive increases of k in 
the ground state by one. From Eq. (fT^ . the average value of the magnetic field, which 
corresponds to the jump k k + 1, follows to be 

m = [{2k + p + l){d)- {A,El,)] . (48) 

Due to the paring interaction, the values {Hk) shift by (— Afc^EQ^) towards higher mag- 
netic fields as compared to the case of normal grains. Since the pair breaking energy 
—AkE^i. decreases with increasing k, the aforementioned shift decreases, too. While in 
normal grains the difference {Hk+i — Hk) between adjacent values (Hk) takes a constant 
value (d) / fi^ for all k, in superconductor grains this difference increases with increasing k 
from (2(d) - {AlEl^)]^^^^) /(2/iB) for A; = to {d)/fXB for k > n^°°\ This means that for 
H > (Hq) the density of the values {H^) on the magnetic-field axis decreases with rising 
H. Consequently, for H > (Hq) the magnetic susceptibility of an assembly of supercon- 
ductor grains should decrease with H, approaching the value x = 1 when H ^ oo. Since 
at magnetic fields H, which are considerably lower than (Hq), the magnetic susceptibility 
of superconductor grains vanishes, we conclude that the afore-described behavior of {H^) 
should result in the appearance of a maximum of x{H, 0) at if ~ (Hq). 

Second, it is worth noticing that the expression — Afc-EofcL=o~'~^-f^+i"~^^' '^^ich determines 
the value Hq for an even grain, coincides with the definition of the spectroscopic gap (cp. 
Ref. Q). As discussed in Ref. the pair-breaking energy — Afc^'g^l^^Q is a decreasing 
function of the spacing Ek+i — £k- Therefore, in superconductor grains, variations of the 
magnetic field Hq with varying the spacing Ek+i — Ek are smaller than those in normal 
grains. In other words, due to the pairing interaction, the values Hq, related to various 
grains, tend to shrink towards {Hq) (for Hk with 1 < k < n^^^\ such a trend is much less 
pronounced, because with increasing k the pair-breaking energy decreases, while variations 
of the spacing Ex+k+p+i — ^K-k increase). The described effect constitutes the second reason 
for the formation of the maximum of x{H, 0) at if ~ (Hq)- 
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The third reason is related to the grains, where at a certain magnetic field Hk,Ak the 
number of doubly occupied levels Ei in the ground state abruptly reduces by a value Ak 
larger than one. At T = 0, each transition of this kind results in an abrupt jump of the 
magnetization by a relatively large value 2/iBAA;. In view of Eq. (jl8j) . it is obvious that 
those transitions lead to an increase x{H, 0) a.t H < [and, correspondingly, in a decrease 
of x{H, 0) at H > Hk] as compared to the case when k changes by one. Since the very 
possibility of these transitions is determined by the pairing interaction, their probability 
generally takes maximum values at /c = 0, when the pairing interaction is the most efficient. 
Therefore, the contributions of those transitions to 0) take maximum values for H < Hq 
enhancing the magnetic susceptibility at if ~ (Hq). 

With increasing k, the number of blocked energy levels in the interaction band increases, 
resulting in a weakening of pairing correlations. At relatively high magnetic fields, when 
the number k becomes larger than n^'^^\ the difference between magnetic susceptibilities 
of superconductor and normal grains disappears. An increase of temperature also tends 
to reduce the above difference. Nevertheless, even for k-^T ~ (d), this difference is still 
observable: while for normal grains the magnetic susceptibility x nionotonously decreases 
with increasing H, for superconductor grains it first increases and then either saturates or 
reaches a maximum and further decreases towards the saturation value x = 1- 

In Figs. El^a) and Efb), the normalized magnetic susceptibility x{H,T) is shown for, 
respectively, normal and superconductor grains with P = 1. Regardless of the value of A, 
the magnetization of an odd grain contains a contribution due to the electron, which is not 
paired in the ground state at H = 0. In the case of T = 0, the spin of this electron is 
fully polarized at arbitrarily weak magnetic fields, resulting in a (5-like peak of the magnetic 
susceptibility at H = 0. Apart from this peak, the behavior of the magnetic susceptibility 
for an ensemble of odd normal grains at T = is described by the equation 



which can be easily obtained using Eqs. ()45|) and PU)) with P = 1. At low magnetic fields, 
the behavior of (x) is mainly determined by the function wi{2^-bH) in Eq. (j^6|) . As follows 
from Eq. the function wf°''''°^(2/iBii) goes to zero when decreasing H. 

For superconductor grains, the behavior of x{H-, T) is qualitatively similar in the cases of 
even an odd grains, except for the aforementioned peak of x{H,T) related to the presence 




(49) 



20 



of an unpaired electron in the ground state of an odd grain at H = 0. Due to the presence 
of this unpaired electron, one of the energy levels Si in the interaction band of an odd grain 
is always blocked. Therefore, in odd grains, the interaction energies are, on average, lower 
and effects, related to pairing correlations, are less pronounced as compared to those in even 
grains. As a result, the maximum of x{H, 0) at non-zero magnetic fields is somewhat lower 
for odd grains than for even grains. One can also see that, as compared to the case of even 
grains, for odd grains this maximum is shifted towards higher magnetic fields. The reason 
for this shift can be easily understood by inspecting the expression on r.h.s. of Eq. ()16|) in 
the cases of P = and P = 1. 

B. Effect of the bare-level statistics on the magnetic susceptibility 

Since an experimental realization of grain assemblies, which would contain only even (odd) 
grains, looks problematic, the case of grain assemblies with equally probable values P = 
and P = 1 seems to be important. Figure ^ displays the plots of x{H,T) for assemblies, 
where even and odd parallelepiped-shaped grains are present in equal proportion. For normal 
grains, the magnetic susceptibility is seen to be a monotonously decreasing function of H. 
For superconductor grains, the magnetic susceptibility as a function of H has two well- 
distinguishable maxima, separated by a deep minimum, at low temperatures {ksT ^ (d)) 
and monotonously increases with magnetic field at ~ (d). 

In order to examine how sensitive are the obtained results to the single-electron energy- 
level statistics, let us consider now the case when energy levels Ei in a grain follow the 
Gaussian orthogonal ensemble (GOE) distribution with the mean nearest-neighbor level 
spacing d. In order to make the results commensurable with those obtained above, we 
assume that the value d varies within the assembly of grains in the same way as in the case 
of parallelepiped-shaped grains. The calculated plots of x{H, T) are shown for normal and 
superconductor grains in Figs.Efa) andEfb), respectively. For a GOE ensemble, probabilities 
of different energy spectra {Si} may differ from each other by many orders of magnitude. 
Therefore, suppression of fluctuations, related to random distributions of energy levels 
in individual grains, is not very efficient even at a relatively large number of grains in 
an assembly. Apart from those relatively large fluctuations, the behavior of x{H, T) for 
superconductor grains is qualitatively the same in the cases of the Poissonian and GOE 
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statistics [cp. Fig. Efb) and Fig. Efb)]. At the same time, the maximum of x{H,0) at 
non-zero H is lower in Fig. Efb) than in Fig. Efb). This is because small values of the 
nearest-neighbor spacing ej+i — Sj are unfeasible in GOE (cf. "pseudogap" ) , as distinct 
from the case of the Poisson distribution of this spacing. Therefore, the interaction energy, 
which has been demonstrated to play a determinative role in the formation of the maximum 
of x{H-,^)-, is, on average, smaller in the case of the GOE statistics as compared to the 
case of the Poissonian statistics Q|. Nevertheless, at T ^ the maximum of x{Hi^) at 
non-zero H exceeds 1.5, while the calculations, performed in Ref. |9|for T = and A = 0.28, 
give the height of the maximum as small as 1.02. We relate this discrepancy to the use 
of purely perturbative results for the interaction energy in Ref. 0. It is noteworthy that 
exact calculations, carried out for a single nanograin with equally spaced energy levels Sj 
[d = 6.6 meV) at temperatures as high as T ~ O.Srf/fce, show a maximum of x{H, 0), which 
significantly exceeds 1.02d/ (2/i|) (see Appendix B). 

As mentioned above, for energy levels described by the GOE statistics, small values 
of interlevel spacing have low probabilities. In this sense, the behaviour of the function 
which describes the probability density for the nearest-neighbor spacing s in a 
GOE ensemble [10], is similar to that of the function That is why the plots of 

x{H, T), shown in Fig. El and calculated for an assembly of both odd and even grains, appear 
resembling those for odd grains with the Poissonian statistics of single-electron energy levels 
(cp. Fig. El and Fig. IH}. Due to low probabilities of small interlevel spacings in normal 
grains with the GOE statistics of single-electron levels, a pronounced minimum appears in 
the low-temperature magnetic susceptibility x as a function of magnetic field at H close to 
zero. In view of this, it is rather the presence of a minimum in the temperature dependence 
of the zero-field magnetic susceptibility 0, 0] and of the maximum in x{H) at i7 ~ (Ho) 
than the presence of a minimum in the low-temperature magnetic susceptibility xiH) at low 
H what can serve to distinguish the magnetic response of superconductor nanograins from 
that of normal nanograins. 
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C. Effect of the grain size and the interaction strength on the magnetic suscepti- 
bility 

Figureini where we go back to the case of parallelepiped-shaped grains, shows the magnetic 
susceptibility x{H,T) for an assembly of superconductor grains (both even and odd) with 
Ko = 1000. When comparing Fig. IHlwith Fig.01^b), which is for grains with two times larger 
(on average) volume, one can see that with reducing the grain size the peak of x{H, 0) at 
nonzero H becomes lower and shifts towards smaller values of fiBH/{d). Of course, the 
absolute value of H, which corresponds to this peak, increases with decreasing the grain size 
(or, equivalently, with increasing (d)). The effect of decreasing the grain size on x{H, T) can 
be qualitatively understood when taking into consideration the following two circumstances. 
On the one hand, a decrease of the grain size leads to an increase of the spectroscopic gap 
[isj l . On the other hand, although the pair-breaking energy rises with reducing the grain size 
(or, equivalently, with increasing (rf)), this rise is slower than an increase of {d). Therefore, 
the relative contribution of the pair-breaking energy to the spectroscopic gap decreases when 
grains become smaller. 

A comparison of Fig. [71 where we plot x(-ff, T) for A = 0.15, with Fig. Efb), which cor- 
responds to the case of A = 0.22, gives an idea of how the magnetic susceptibility x{H,T) 
of superconductor grains changes with weakening the pairing interaction. When consider- 
ing the magnetic susceptibility as a function of the dimensionless variables ^bH/ (d) and 
k-BT/{d), the effect of decreasing A on x looks very similar to that of decreasing the grain 
size [cp. Figs. Hfb), El andHj. This similarity is natural, because both a decrease of A 
and a decrease of the grain size reduce the ratio of the pair-breaking energy to the average 
nearest-neighbor level spacing (d). 

It is worth noting that the behavior of x{H,T), described throughout this section, is 
sufficiently stable with respect to an increase of the grain-size dispersion 6 within an assembly 
of grains. As mentioned above, due to an extremely high sensitivity of the energies Si to the 
grain shape and size, even within a set of grains with a relatively weak dispersion of sizes 
{6 = 0.025), one finds a vast variety of different single-electron energy spectra. Therefore, a 
moderate increase of 6 almost does not affect the dispersion of pairing characteristics over an 
assembly of grains. For the magnetic susceptibility of grains with 6 = 0.2, our calculations 
give results, which practically coincide with those in the case of 5 = 0.025 (just for the 
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reason of such a close coincidence, we do not show the plots for S — 0.2). 

IV. CONCLUSIONS 

We have shown that a simple approximation developed here for the interaction energy in 
superconductor grains is applicable over a wide range of the grain sizes and of the dimension- 
less interaction parameter A values at arbitrary distributions of single-electron energy levels 
in the interaction band. This approximation allows us to analyze the magnetic response 
of relatively large assemblies of nanograins even at relatively high temperatures and mag- 
netic fields when the thermal energy and the magnetic energy /ibH are larger than the 
mean interlevel spacing in grains, d, so that the number of many-electron states, involved 
in calculations, can exceed 10^ per grain. 

The overall profiles of the magnetic susceptibility of superconductor nanograins as a 
function of magnetic field and temperature are demonstrated to be qualitatively different 
from those for normal grains. At T — > 0, the normahzed magnetic susceptibility x{H) of 
superconductor grains is characterized by a well pronounced maximum at a certain value 
of 2fiBH, which is close to the width of the spectroscopic gap. Such a maximum or its 
reminiscence in the form of an increase of x with H at low magnetic fields persist even 
at k-sT ~ d, both for even and odd grains as well as for their mixtures. The analyzed 
distinctions between the magnetic response of superconductor grains and that of normal 
grains are shown to be qualitatively the same both for the GOE statistics of single-electron 
energy levels in grains and the Poissonian one. This "universality" seems important when 
taking into consideration the fact that an appropriate choice of statistics to describe the 
energy-level distribution in concrete nanostructures may be not self-evident. The aforemen- 
tioned features are sufficiently stable also with respect to variations of the average value of 
the grain size and its dispersion over an assembly of nanograins. These features are shown 
to be quite pronounced for the ultra-small superconductor grains, where both the spacing 
d and its variations within an assembly of grains significantly exceed the superconducting 
gap A, so that pairing correlations can hardly be detected through the tunneling spectra. 
Our results imply that measuring the magnetic-susceptibility profiles, x{H,T), can provide 
a sensitive probe of the pairing interaction in those grains. 
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APPENDIX A: EXAMPLE: A GRAIN WITH DEGENERATE SINGLE- 
ELECTRON STATES IN THE INTERACTION BAND AT H = 

Here we analyze the magnetic response of a grain, where in the absence of a magnetic field 
all the single-electron states in the interaction band have one and same energy: Si = Ek+p 
for 2 G X. Assuming that the dimensionless interaction parameter is small, A <^ 1, we restrict 
magnetic fields and temperatures under consideration by the inequalities (ibH <^ hwj) and 
ksT <^ huj-[). Under these conditions, excitations, which are related to a change of occupation 
numbers for the single-electron states outside the interaction band, can be neglected when 
calculating the magnetization of the grain. Therefore, here we consider only many-electron 
states, where all the energy levels Ei below the interaction band are doubly occupied, while 
all the levels above the interaction band are empty. 

Using Eqs. Q and ()33|). the energy of 2n*^°°^ +P electrons, which reside in the interaction 
band, can be written as 

E,,, = (2n(°°) + P) Ek^p - - k - g) {l - P - - k - g + l) 

-{2k + P -2l)fiBH. (Al) 

The degree of degeneracy of the energy level Egki is described by the expression 



where the last factor gives 22 1 the number of (n^^'^-' — A;)-pair states with a given value of g. 
For a state with given k and /, the spin magnetization is described by Eq. (fTUI) . Inserting 
Eqs. |T(H) . (jAlll and ()A2j) into Eq. |TT| and performing summation over Z, we obtain for the 
magnetization of the grain: 



\ / ^max \ ^max 

-t)[Y.^A $^(2A; + P)/^fc, (A3) 



, A:=0 / A;=0 
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where 




exp 



2k{I -P~k + l) 



P-2k 



It 

) exp 



2g{I-P-2k-g + l] 
It 



(A4) 



(/_P_ |j_P_2r2(oo)|) 



(A5) 



h = fisH/i^^D), and t = k-QT/{\f}WD). 

As discussed in Subsection II. B, the magnetization at T = is determined by the polar- 
ization of electron spins in the ground state, which must be found by minimizing the energy 
i?ofco as a function of k. From Eq. (jAlj) one finds A^ii^ofco = —^\fujjY)/I < 0. This means 
that the ground-state energy is given by an edge minimum of -Eofco at /c = or at A; = /Cmax- 
From the equation i^ooo = -^ofcmaxO we find that the transition from the "maximal-pairing" 
ground state (with = 0, n^'^^^ pairs, and magnetization Mqo = I^bP) to the "maximal- 
magnetization" ground state [with k = fcmaxj {n^'^^^ — A;max) pairs, and the magnetization 
fcmaxO = /^B(2/i;max + P)] occurs at the magnetic field 

XhUY) {I - P - fcmax + 1) 



H0,kn 



I^bI 

Xhwiy {I-P+\I -P- 2n(°o) I + 2) 
2iI^I 



(A6) 



The field -ffo.fcmax is determined by the pair-breaking energy —A^Elf^ averaged over fcmax 
pairs broken in the course of the above transition. Since for the system under consideration 
the second forward difference A^_Eofco/ coincides with — A^i^Qj^., the inequality AlEoko < 0, 
derived above, implies that the energy —AkE^k decreases with increasing k. Therefore, the 
value -ffo,fcmax decreases with increasing fcmax or, equivalently, with reducing deviations from 
half filling, which are described by the quantity | / — P — 2rS^'^'> | . 

The field ifo,fcmax can be considered as the zero-temperature critical magnetic field for the 
phase transition between the superconducting and normal states of the grain. At T = 
this transition is characterized by a 5-like peak of the magnetic susceptibility as a function 
of H . At non-zero temperatures, the aforementioned peak is smoothed out. For T 7^ 0, we 
define the "critical magnetic field" as the value of iJ, which corresponds to the maximum 
of the magnetic susceptibility at a given T. The plot of the "critical magnetic field" versus 
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temperature can be considered as a "phase boundary" between the superconducting and 
normal states of the grain. An example of those "phase boundaries", calculated with the 
use of Eq. flA3|) for even grains with J = 50 and different filling, is shown in Fig. |S1 While 
the critical magnetic field at T = increases with increasing deviations from half filling, the 
"critical temperature" at if —> is seen to demonstrate the opposite trend. This is because 
the "critical temperature" at if — is related to the total energy of breaking /cmax electron 
pairs rather than to the average breaking energy per pair, which determines iicfc^ax- 

For odd grains with half filling of the interaction band, an analysis of the magnetization 
at magnetic fields shghtly above iio,fcmax low temperatures {kBT <^ XHub) reveals a well- 
pronounced re-entrant behavior [see Fig.|Ul^a)]. With increasing temperature, magnetization, 
which equals Mmax at T = 0, rapidly reduces downwards the value Mqo, characteristic for the 
"maximal pairing" state. This re-entrant behavior is due to the fact that at /cbT ^ 
the magnetization in the vicinity of iio,fcmax is mainly determined by the interplay between 
the occupation probabilities for the "maximal magnetization" and "maximal pairing" states. 
As seen from Eq. ()A2|1 in the case of P = 1 and 2n*^'^'^^ = I — P, the energy level i^ofcmaxO 
is non-degenerate, while the level i?ooo is multiply degenerate (Jooo = I)- Therefore, for 
H ^ iio,A:max ^iid k^T ~ |i?ofcniaxO ~ -^000 1 ) thcrmodyuamic average magnetization is close 
to Moo even at magnetic fields above H^^k^^^- As seen from Fig.Efb), such a re-entrance does 
not occur in even grains with half filling, where both energy levels f^ofcmaxO aiid i^ooo are non- 
degenerate. When deviating from half filling, the ratio JoA;maxo/'^oo becomes significantly 
larger than 1. As illustrated by Fig. Efc), in those distinction of the temperature 

from T = results in a rapid increase of magnetization at ii < iio.fe^ax from Moo towards 

When increasing i, the ranges of H and T, which correspond to the afore-described re- 
entrant behavior in odd grains with half filling, gradually reduce (see Fig. [TUI) . The reason 
for such a reduction is as follows. As we have seen above, the re-entrant behavior is related 
to a dominant role of the energy levels i?ooo and i^ofcmaxO in determining the thermodynamic 
average magnetization M. However, an increase of i results in a rapid increase of degeneracy 
Jgki for higher excited states and hence in an increasing contribution of those states to M 
at T 7^ 0. Consequently, the relative contribution of the "maximal pairing" states to M 
decreases. The larger is i, the narrower are the range of magnetic fields above iio.fcmax and 
the temperature range, where the "maximal pairing" states play a determinative role for the 
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thermodynamic average magnetization. 



APPENDIX B: ON THE ACCURACY OF THE PROPOSED APPROACH 

In order to assess the accuracy of our approximate approach, let us compare the results, 
obtained within this approach, with those of numerically exact calculations. The latter 
calculations are time-consuming and can be actually performed only for a relatively small 
number of relatively small grains. In Fig. ^2 we compare the approximate and numerically 
exact results for the zero-temperature magnetic susceptibility of superconductor grains with 
Ko = 1000. Even for those very small grains the computation time for the exact calcula- 
tions exceeds that for the corresponding approximate calculations by more than 3 orders of 
magnitude (the above value dramatically increases with rising grain sizes). Therefore, the 
number of grains in the exact calculations is limited to the values, which are not yet large 
enough for a sufficient suppression of fluctuations, related to seemingly random distribu- 
tions of the energy levels {si} in individual grains. Nevertheless, we can see from Fig. 10 
a reasonable agreement between the approximate and exact results. Since our approximate 
approach somewhat overestimates the interaction energy [15], the maximum of the approx- 
imate curve x{H) is higher and lies at stronger H as compared to the exact result. The 
approximate curve shows also a splitting of this maximum. At T — > 0, a similar splitting is 
seen also in Figs. infb),ini andUl From the analysis of the approximate results for assemblies 
of grains with a definite parity of the number of electrons, we may draw a supposition that 
this splitting is related to the parity effect on the position of the maximum of x{H) (namely, 
for odd grains this maximum is shifted toward higher magnetic fields as compared to the 
case of even grains). However, relatively large fiuctuations of the numerically exact curve 
do not allow us to definitely conclude whether such a splitting should be actually observable 
or its appearance/magnitude is an artifact of the approximate approach. 

Numerically exact calculations of the magnetic susceptibility at non-zero temperatures 
are practically possible only for single grains. As an example, we consider an even grain 
with equidistant energy levels ej spaced hj d = 6.6 meV. At k^T ~ {d), calculations of 
x{H, T) for such a grain involve about 5 x 10^ many-electron states (this value gives an idea 
of difficulties in extending the exact calculations to an assembly of grains). In Fig. IT^ a). 
the results for a normal grain (A = 0) are shown. Positions of peaks of x{Hi 0) correspond 
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to those given by Eq. (fT^ for the case of the equidistant energy levels ej: 

ft = i?^, (Bl) 

Figures IT^ b) and lT^ c) display, respectively, the approximate and exact results for A = 0.22. 
Like in the above case of T = 0, we see a satisfactory agreement between these results. As 
clear from Figs. IT^ b) and IT^ c). in the grain under consideration, the number of doubly 
occupied energy levels Ej in the ground state, K — k, decreases consecutively by one with 
increasing H. In this grain, the pair-breaking energy {—AkEli^) and, correspondingly, the 
difference between the field H^. for the superconductor grain and that for the normal grain 
are significantly larger for k = than for k > 1. As a result, for the superconductor grain 
the distance Hi — Hq between the first two peaks of x{H, 0) appears significantly smaller 
than the distances Hk+i — Hk for the further peaks (with k>l). This effect is a single-grain 
analog of the behavior of x{Hi 0) in assemblies of superconductor grains, which is discussed 
in Subsection III. A: due to the pairing interaction, the magnetic susceptibility of an even 
superconductor grain almost vanishes at low magnetic fields and takes the maximum values 
when the magnetic energy 2ii^H becomes close to the spectroscopic gap. 
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FIG. 1: Condensation energy, calculated for parallelepiped-shaped grains with the aspect ratio 
Ix '■ ly '■ h = ^ '■ + v) ■ ~^ 27/), as a function of the parameter rj. The volume of grains is 
fixed to keep a constant number of electrons in a grain, 2K = 4000. The dimensionless interaction 
strength, A, is taken to be 0.22. The dotted curve represents the numerically exact results Q]. 
The solid curve displays the results based on the use of Eq. (|l0|l . 
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FIG. 2: Normalized magnetic susceptibility as a function of temperature and magnetic field for two 
assemblies of parallelepiped-shaped even grains: (a) normal grains (A = 0) and (b) superconductor 
grains (A = 0.22). For each of these assemblies, the total number of grains is 6 x 10^, the parameter 
Kq is 2000, and the average nearest-neighbor level spacing in the interaction band is (d) = 3.1 meV. 
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FIG. 3: Normalized magnetic susceptibility as a function of temperature and magnetic field for two 
assemblies of parallelepiped-shaped odd grains: (a) normal grains (A = 0) and (b) superconductor 
grains (A = 0.22). For each of these assemblies, the total number of grains is 6 x 10^, the parameter 
Kq is 2000, and the average nearest-neighbor level spacing in the interaction band is (d) = 3.1 meV. 
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FIG. 4: Normalized magnetic susceptibility as a function of temperature and magnetic field for 
two assemblies of parallelepiped-shaped grains: (a) normal grains (A = 0) and (b) superconductor 
grains (A = 0.22). For each of these assemblies, the total number of grains is 1.2 x 10^, even and odd 
grains are mixed in equal proportion, the parameter Kg is 2000, and the average nearest-neighbor 
level spacing in the interaction band is (d) =3.1 meV. 



34 




FIG. 5: Normalized magnetic susceptibility as a function of temperature and magnetic field for 
two assemblies of grains with the GOE statistics of single-electron energy levels: (a) normal grains 
(A = 0) and (b) superconductor grains (A = 0.22). For each of these assemblies, the total number 
of grains is 1.4 x 10^, even and odd grains are mixed in equal proportion, and the average nearest- 
neighbor level spacing in the interaction band is (d) = 3.1 meV. 
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FIG. 6: Normalized magnetic susceptibility as a function of temperature and magnetic field for 
an assembly of parallelepiped-shaped superconductor grains with A = 0.22. The total number of 
grains is 1.2 x 10^, even and odd grains are mixed in equal proportion, the parameter Kq is 1000, 
and the average nearest-neighbor level spacing in the interaction band is (d) = 6.6 meV. 
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FIG. 7: Normalized magnetic susceptibility as a function of temperature and magnetic field for 
an assembly of parallelepiped-shaped superconductor grains with A = 0.15. The total number of 
grains is 1.2 x 10^, even and odd grains are mixed in equal proportion, the parameter Kq is 2000, 
and the average nearest-neighbor level spacing in the interaction band is {d) = 3.1 meV. 
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FIG. 8: "Phase boundaries" for even superconductor grains with 2n^^^^ electrons distributed over 
27 single-electron states in the interaction band, all these states possessing the same energy in the 
absence of a magnetic field. 
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FIG. 9: Magnetization of superconductor grains, where 21 single-electron states in the interaction 
possess the same energy at = 0, as a function of temperature and magnetic field. The number of 
electrons in the interaction band is 2n^^'^') + P. Heavy dashed lines show the corresponding "phase 
boundaries". Panel (a): an odd grain with half filling of the interaction band. Panel (b): an even 
grain with half filling of the interaction band. Panel (c): an even grain with a slight deviation from 
half filling of the interaction band. 
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FIG. 10: "Phase boundaries" for odd superconductor grains with different number 21 of single- 
electron states in the half-filled interaction band, all these states possessing the same energy in the 
absence of a magnetic field. 
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FIG. 11: Normalized magnetic susceptibility as a function of magnetic field for parallelepiped- 
shaped superconductor grains with A = 0.22 at T = 0. The dashed curve shows the numerically 
exact results for an assembly of 4 x 10^ grains. The solid curve displays the results, based on the 
use of Eq. (|4(H) . for an assembly of 6 x 10^ grains. In both cases, even and odd grains are mixed in 
equal proportion, the parameter Kq is 1000, and the average nearest-neighbor level spacing in the 
interaction band is {d) = 6.6 meV. 
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FIG. 12: Normalized magnetic susceptibility as a function of temperature and magnetic field for 
single even grains with the equidistant single-electron energy levels spaced hy d = 6.6 meV. Panel 
(a) shows the results for a normal grain (A = 0). For a superconductor grain with A = 0.22, the 
approximate results, based on the use of Eq. (|4U|) . and the numerically exact results are displayed 
on panels (b) and (c), respectively. 
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